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Abstract 

It has been recently claimed that the symmetry group S4 yields to the Tri-bimaximal 
neutrino mixing in a "natural" way from the group theory point of view. Approving 
of this feature as an indication, we build a supersymmetric model of lepton and quark 
masses based on this family symmetry group. In the lepton sector, a correct mass 
hierarchy among the charged leptons is achieved together to a neutrino mass matrix 
which can be diagonalized by the TVi-bimaximal pattern. Our model results to be 
phenomenologically unequivalent with respect to other proposals based on different 
flavour groups but still predicting the Tri-bimaximal mixing. In the quark sector a 
realistic pattern for masses and mixing angles is obtained. The flavour structures 
of the mass matrices in both the sectors come from the spontaneously symmetry 
breaking of S4, due to several scalar fields, which get non-zero vacuum expectation 
values. A specific vacuum alignment is required and it is shown to be a natural 
results of the minimization of the scalar potential and, moreover, to be stable under 
the corrections from the higher order terms. 
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1 Introduction 



By now there is convincing evidence that the solar and the atmospheric neutrino anomalies 
can be explained by the neutrino oscillations. The Am 2 values and mixing angles are known 
with good accuracy [1-3]. The latest best values for Am 2 are Am^ m ~ 2.4 x 10~ 3 eV 2 
and Am 2 oZ ~ 7.7 x 10~ 5 eV 2 . For the mixing angles, two are large and one is extremely 
small: the atmospheric angle 6*23 is compatible with a maximal value, but the accuracy 
admits relatively large deviations, indeed at 2a errors it is 0.366 < sin 2 #23 < 0.602 with 
central value 0.466; the solar angle #12 is large, 0.278 < sin 2 9 12 < 0.352 with central value 
0.312, but about 5a errors far from the maximal value; the reactor angle #13 is strongly 
bounded and at present it has an upper limit of sin 2 #13 < 0.036. We underline that there 
are contrasting indications for a vanishing value of the reactor angle: in [2] there is a 
suggestion for a positive value which, at 1.6a, is sin 2 9 13 ~ 0.016 ± 0.010, while in [3] the 
authors find a best fit value consistent with zero within less than la. Therefore we need 
a confirmation by the future experiments like DOUBLE CHOOZ [4], Daya Bay [5] and 
MINOS [6] in the v e appearance channel. 

From the theoretical point of view, the developments about neutrino masses and mixing 
angles cannot satisfy: there is a so large number of existing models, that can be interpreted 
as a lack of a unique and compelling theoretical picture. However a series of models 
based on some discrete non-Abelian groups seems to be extremely attractive due to their 
predictions: indeed it is possible to achieve as the lepton mixing matrix the Tri-bimaximal 
(TB) pattern [7], 
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which represents a very good approximation of the experimental data [2], providing the 
following values for the mixing angles: 

sin 2 6™ = sin 2 9™ = 1/2 sin 2 9™ = 1/3 . (2) 

It is just the TB pattern which can suggest the type of symmetry that best describes the 
lepton mixings: it is a very well known result [8] that a maximal value for the atmospheric 
angle can be recovered only with a non exact symmetry; explaining the indication for a non 
vanishing, but still very small, value for # 13 , it is necessary to provide the TB pattern at 
the leading order (LO), invoking corrections from the higher order terms; the solar angle is 
predicted to be very close, less than 2°, to the measured value and therefore the corrections 
has to be relatively small. As a result, a realistic lepton flavour symmetry has to be broken 
at a certain level, predicting at the LO the TB pattern and providing corrections at the 
next-to-the-leading-order (NLO) at most of about 9 2 c ~ 2°, where 9 C stands for the Cabibbo 
angle, which is a convenient hierarchical parameter for both the sectors. 

There is a series of models based on the symmetry group A4 [8-15], which are extremely 
attractive from this point of view, fulfilling all the previous requirements. A 4 is the group 
of the even permutations of four objects and has 12 elements and four irreducible repre- 
sentations, which are three singlets, 1, V and 1", and one triplet 3. These models manage 
in deriving the TB mixing by assuming that the A4 symmetry is realized at a very high 
energy scale A and that leptons transform in a non trivial way under this symmetry. After- 
ward the group is spontaneously broken by a set of scalar multiplets 0, the flavons, whose 



vacuum expectation values (VEV) receive a specific alignment. It is a non trivial task to 
explain how to get the expected vacuum alignment in a natural way and we consider it a 
fundamental requirement for a competitive model. Moreover the TB mixing is corrected by 
the higher order terms by quantities of the order of (0)/A < 1 and as a result the reactor 
angle is no longer vanishing and becomes proportional to (0)/A. 

The common aspect of many of these projects is the structure of the neutrino mass 
matrix. The most general mass matrix for the neutrinos which can be diagonalized by the 
TB mixing is the following 



and it is /i t invariant, yielding to a maximal atmospheric mixing angle and to a 
vanishing 6^3 , and it satisfies the relation m Ull + m Ul3 = m U22 + m U23 , which gives the Tri- 
maximal solar angle [17]. For # 13 = there is no CP violation from the Dirac phase, and 
there are only Majorana phases. If we disregard them, we can restrict our considerations 
to real parameters. Usually this pattern can be obtained constructing the lagrangian in 
such a way that the usual Weinberg operator, which we can write as ££ implying £h u £h u , is 
forbidden at the leading order, but appears only at higher orders with additional flavons. 
Most of the models based on the flavour symmetry, are characterized by b = Q; 
however, different realizations with other relations between a, b and c have been studied, 
see for instance [15]. In the pattern with 6 = 0, the factors a in eq.([3]) come from the term 
££F\ and the factors c from ££F 3 , where Fi and F 3 are flavons transforming respectively as 
a singlet 1 and as a triplet 3 of v4 4 . Presenting the same flavour structure for the neutrino 
mass matrix, it is extremely difficult to distinguish one model from all the others by the 
use of only observables connected to the neutrino oscillations. Some improvements in this 
direction has been recently performed in [18], where the authors develop an analysis on 
some lepton flavour violating processes, which can be tested in the future experiments. 

Moreover the great difficulty of this kind of models is to describe correctly the quark 
sector. First of all the quark mixing matrix is completely different from its lepton coun- 
terpart: the first shows little angles and, in the contrary, the second presents two large 
angles. As a result, while the lepton mixing matrix can be fairly achieved through a dis- 
crete flavour symmetry, the quark mixings seem to be better described by some continuous 
symmetry, like U(2) [19]. Indeed, according to the left and right-handed quark representa- 
tion assignments, a discrete non-Abelian flavour symmetry tends to predict no mixing at 
all in the quark sector, Vckm = 1, or too large mixing angles. On the other hand, the re- 
sults obtained by the U (2)-based models for the quarks suggest that the use of the doublet 
representation in the quark sector should help in describing quark mixing. However, this 
possibility is prevented in the y4 4 -based models, since there are not doublet representations. 

lr rhe pattern with 6 = can be obtained with A4 as flavour symmetry only with a particular flavon 
spectrum: the singlets 1' and 1" of A4 have not to couple to the term It. In fact with the left-handed 
leptons, £, transforming as triplet of A4, we can switch on the entries corresponding to b by coupling 
II to the flavons Fv and singlets 1' and 1" respectively, as can be checked by looking in [10] and 
already underlined in the Appendix of [16]. Moreover, in this case, the two couplings give two distinct 
contributions and as a result the mass matrix is not diagonalizable by the TB pattern any more. To be 
diagonalized by this mixing scheme, it is necessary to impose the condition y\(Fy) = 2/2 (-Fi" ), where yi 
are the coupling constants of the two operators. However, keeping the model as natural as possible, i.e. 
without fine-tuning, it is necessary to prevent the couplings of the two singlets Fy and Fy with ££. The 
commonly used solution consists in not introducing such flavons. 




(3) 



The solutions which have been proposed consist in the possibility of add several Z n sym- 
metries [12], in order to suppress the unwanted terms, or in adopting a larger group, which 
manages in reproducing the structure of in the lepton sector and possesses some doublet 
representations useful to describe quarks, like for example the discrete group T' [20,21]. 
In our opinion, a good candidate to be the flavour symmetry group describing leptons and 
quarks has to be as small as possible and has not to need of numerous additional elements 
in order to reproduce correct fermion masses and mixing angles. Following this prejudice, 
we looked for a proposal which manages in describing both the sectors in a realistic way, 
keeping as simple as possible the symmetry content. 

It has been recently claimed [22] , through group theoretical arguments, that the minimal 
flavour symmetry naturally related to the TB mixing is sj^l [23,24]. The group S4 is the 
group of the permutations of four objects and it has 24 elements divided into five irreducible 
representations: two singlets li and I2, one doublet 2 and two triplets 3i and 32 (a more 
detailed description of 5*4 can be found in Appendix A). We approve of the result in [22] 
as an indication and we present a model based on the discrete non-Abelian symmetry 
group £4, which predicts the TB mixing in the lepton sector and a CKM matrix close 
to the experimental one in the quark sector (the group S4 has already been studied in 
literature [25], but with different aims and different results). Moreover we introduce an 
additional Z 5 symmetry, which plays a similar role of the total lepton number avoiding 
some dangerous terms, and a continuous U(1)fn [26], that helps to provide the correct 
fermion hierarchies. S4 contains as a subgroup A4, but it has a doublet representation, 
which can be used in order to describe quarks. It has the same number of elements of T", 
but the representations are different: in particular T' can derive only the same neutrino 
mass matrix of the ^U-based models. On the other hand, the mass matrix which can 
be constructed in a S^-based model is exactly that one in eq.([3]) and therefore it is more 
general with respect to the previous case. From this point of view we can say that T' 
constraints the neutrino sector in a stricter way than 5*4. We underline that the model 
cannot be embedded into a GUT context, because of the different transformation properties 
of leptons with respect to quarks H. 

The presence of the doublet representation, not only represents the new expedient in 
order to describe the quark sector, but also introduces a new feature in the neutrino mass 
matrix: indeed the terms which contribute to m u are MFi, i£F 3 and the new i£F 2 , where 
F2 represents a flavon transforming as a doublet 2. In eq.®, this last contribution is 
represented by the term b. This result corresponds to the neutrino mass matrix in [24]: 
however the presence of three parameters in order to describe three masses prevents any 
predictions on the neutrino hierarchy type. For these reasons we conclude that the S4- 
based model in which a singlet F±, a doublet F 2 and also a triplet F 3 couple to 11 is 
not phenomenologically interesting. However it is not restrictive to construct a model in 
which only a singlet and a doublet contribute to the neutrino mass matrix, but in this 
case mi = 7773 and it would be spoiled out by the experimental observations. Moreover 
it is possible to think about a model in which only a singlet and a triplet contribute to 

2 We agree with the conclusions of the group theoretical analysis, but, in our opinion, it must not 
be considered a constraint for the model realization: from the model building point of view, the most 
economical realization which naturally provides the TB pattern as the neutrino mixing matrix is based on 
the A4 symmetry group. 

3 When we were completing our work, the following paper appeared [23] , in which the authors present 
a model based on the symmetry group SU(5) x S4. However in this model it is not possible to explain 
completely the VEV alignment and as a result the mass hierarchies and some mixings have to be fine-tuned. 



the neutrino mass matrix: we have verified that such a model can be built, with a natural 
vacuum alignment. This model provides exactly the neutrino mass matrix with 6 = and 
therefore it has the same predictions in the lepton sector as of the y^-based models. For 
this reason in this paper we study the case in which only a doublet and a triplet couple to 
the term tl and as a result we get an unusual neutrino mass matrix 

(2c b — c b — c \ 
b-c b + 2c -c \ (4) 
b-c -c b + 2c J 

which can still be diagonalized by the TB mixing. This new pattern provides different 
predictions for the Oz/2/3-decay and thus this model can be distinguished from all the oth- 
ers which predict the TB mixing, just looking at some observables related to the neutrino 
oscillations. 

In the following we first provide a phenomeno logical analysis of the new neutrino mass 
matrix, underlining the connections with the Oz/2/3-decay Subsequently, in section 3, we 
present the model which naturally develops the TB mixing in the lepton sector and an ac- 
ceptable CKM matrix in the quark sector in addition to realistic mass hierarchies between 
all the fermions. In section 4, we show how to get in a natural way the special vacuum 
alignment, used throughout the paper. In section 5, we present a study on the corrections 
introduced by the higher order terms. Finally, we summarize the results in the conclusions. 
Details on the group S4, like the conjugacy table, the complete list of the elements in a 
particular basis of the generators and the respective Clebsch-Gordan coefficients, can be 
found in Appendix A. The complete NLO analysis of the vacuum stability is presented in 
Appendix B. 



2 Phenomenological Analysis 

The neutrino mass matrix in eq.(j4]) can be diagonalized by the TB mixing and the eigen- 
values are given by 

v 2 

m dia 9 = ( 3c _ 6; 26, 3c + 6)^ . (5) 
We can now write the neutrino oscillation parameters Am^ m and Am^ as follows: 

Am 2 atm =| m , 3 | 2 -K i | 2 = 12|6||c|cosC^ (6) 

AmL = KJ 2 " KJ 2 = 3|6| 2 - 9|c| 2 + 6|6||c| cosC^ (7) 

where the angle ( is the relative phase between 6 and c. This phase is related to the 
Majorana CP phase 021, which is defined as follows 

E/^UrB-diag^e^e' 2 * 1 ) . (8) 

We can express |6| and |c| as functions of Am 2 tm , Am 2 ol and ( and as a result we get 
constraints on the type of the neutrino spectrum, on the value of the lightest neutrino 
mass and on the 0z/2/3 parameter \m ee \ directly from the experimental data. In fig.(CE|) 



on the left, we plot | function of the lightest neutrino mass eigenstate, m v \ in 

the normal hierarchy (NH) case and m v z in the inverse hierarchy (IH) one. On the right, 
we present | function of the Majorana phase 021- We observe from this last 

plot that considering the Heidelberg-Moscow [27] experiment, which provides the lowest 
present bound on \m ee \ of about 0.35 eV, the exact CP-conserving Majorana phase 02i = 
is excluded in our model. Moreover, from fig.(CE]) on the left, we conclude that the NH region 
falls in the quasi Degenerate Case (DC) band and therefore we cannot speak properly of 
NH in this model. 
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Figure 1: On the left it is plotted \m ee \ as a function of the lightest neutrino mass, m Vl in red in the case 
of the NH and m„ 3 in blue in the case of the IH. The light colored bands represent the possible regions 
considering only the exact TB pattern, while the dark colored ones are the predictions of our model. The 
present bound from the Heidelberg-Moscow experiment is shown in dark gray and the future sensitivity 
of CUORE (~ 15 meV), Majorana (~ 20 meV), and GERDA (~ 90 meV) experiments are represented 
by the horizontal dashed lines, while the future sensitivity of 0.2 eV of KATRIN experiment is shown by 
the vertical dashed line. On the right, \m ee \ as a function of the physical Majorana phase a2\- in red the 
NH case and in blue the IH one. The dark gray region and the dashed horizontal lines corresponds to the 
previous plot. 



Restricting our discussion to the IH case, we find a lower bound for the 0z/2/3 parameter, 
\m ee \ > 14.4 meV, for the lightest neutrino mass, \m u:i \ > 0.72 meV, and for the sum of 
the neutrino masses, Y2i \ m i\ > 89 A meV, all of them corresponding to £ = 0. Moreover 
we have a prediction for \m ee \ in function of Am^ im , r = Am^/Am^ and ( 

-(1 + r)Am 2 atm + yj (AmL) 2 cos 2 C (3(r - l) 2 + (r + l) 2 cos 2 Q sin 2 C . 

We observe that our model can be distinguished from that one in [8,10], based on A4, 
and that one in [20], based on T', looking to the lower bound on |m ee |: in fact those 
models predict a lower bound for |m ee |, which is about 0.005 eV. Our predictions are 
quite close to the future experimental sensitivity, which are expected to reach the values of 
0.090 eV [28] (GERDA), 0.020 eV [29] (Majorana), 0.050 eV [30] (SuperNEMO), 0.015 eV 
[31] (CUORE) and 0.024 eV [32] (EXO). 
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3 The Model 



The discrete group S4 is given by the permutations of four objects and it is composed by 
24 elements. It can be defined by two generators S and T that satisfy 

S 4 = T 3 = (ST 2 ) 2 = 1. (9) 

The three relations reported above directly indicate which are the discrete Abelian sub- 
groups of S4: Z 4 , Z 3 , Z 2 respectively Indeed the 24 elements of S4 belong to five classes 
reported in Appendix A: the elements of C 2} 4 define two different sets of Z 2 subgroups of 
5*4, corresponding to S 2 and ST 2 respectively, those of the class C4 a set of Z 3 Abelian 
discrete symmetries associated to T and those belonging to C5 a set of Z4 Abelian discrete 
symmetries corresponding to S. From the three relations that define the group 1S4 we 
see that it contains also a non-Abelian subgroup, S3. Indeed defining S' = S 2 and using 
S 2 TS 2 = T 2 we get the relations that define S3, namely 

T 3 = S' 2 = (S'T) 2 = 1 . (10) 

Furthermore, S4 presents 5 irreducible representations: two singlets, li, I2, one doublet, 2, 
and two triplets, 3i and 32- All the technical details are reported in Appendix A. 



3.1 The Lepton Sector 

In this part we illustrate the model in the lepton sector, predicting an exact TB mixing 
at the LO and a realistic charged lepton mass hierarchy, by the use of flavour group Gf 
in addition to the gauge group of the SM. The complete flavour group is Gf = S4 x Z 5 x 
U(1)fn, where the three factors play different roles: the spontaneous breaking of S4 down 
to its subgroup Z 2 x Z 2 in the neutrino sector is directly responsible for the TB mixingQ; 
the Z5 factor plays a similar role of the total lepton number, avoiding some dangerous 
terms, and, together to the U(1)fn, is responsible for the hierarchy among the charged 
fermion masses. In table HJ we can see the lepton sector fields of the model and their 
transformation properties under Gf. We treat the model in a supersymmetric scenario, 
because the minimization of the scalar potential is simplified, but this is not a constraint 
from the construction of the model itself. 

The superpotential for the leptons can be written as 



wi 



i=l 



E {^ e Wh d + %^rj)'h d + \r\t^h d (11] 



where 



w v = |(WM + ^(thJKA) (12) 



X = {W V , iPw, AA£', A<p?} (13) 



4 This breaking is extremely unusual, indeed the common preserved subgroup is Zi- Here Z2 x Z2 
provides the same flavour structure for the neutrino mass matrix as Zi in the A4 based models and it is 
associated to one element of the class C2 and one of the class C4. The complete list of the elements are 
present in the Appendix A. 
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Tabic 1: Transformation properties of the matter fields in the lepton sector and 
of all the flavons of the model. We distinguish the flavon fields on their role and 
thus we can consider tjj and rj mainly connected to the charged lepton sector 
and A and ip to the neutrino sector. All these fields together to £' are present 
in the quark sector. The FN field, 9, provides the correct mass hierarchy. 



using (. . .) to refer to the contraction in li and (. . .)' to the contraction in 1 2 . It is 
interesting to underline that the first contributions containing e c would be 

^e%iAA)'h d + ^e%£Acpyh d , (14) 

which would dominate with respect to the terms in eq. ffTTT) . However an explicit compu- 
tation will show that these two terms are vanishing, once we assume that the flavons get 
this specific VEV: 



1 




(A) = 1 Ua {<p) = ( x ) v v 



1 



(15) 



(O = v e (9) = v e 

We will demonstrate that this particular VEV alignment is a natural solution of the scalar 
potential in the following sections; moreover we will see that all the VEVs are of the same 
order of magnitude and for this reason we will parameterize the ratio VEV/ A by the 
parameter u. The only VEV which originates with a different mechanism with respect to 
the others is Vg and we indicate the ratio vg/A by the parameter t. 

With this setting, in the basis of canonical kinetic termg^l, the mass matrix for the charged 
leptons is (rri£ ~ R C L, m u ~ L T L) 

{1) uH yi 2) u 2 t yi 3) u 2 t 




m e = | y^u | uv d (16) 
y T 

where the yi are the result of all the different contributions of the y e ^. For the neutrinos 



5 It has been shown in a series of papers [33] that the corrections, from the transformations needed to 
move in the basis of canonical kinetic terms, appear at most as NLO deviations. 



we get the following mass matrix, which is exactly diagonalized by the TB pattern, 



2c b — c b — c 
i — c b + 2c — c 
b — c —c b + 2c 



m u = | b-c 6 + 2c -c I (17) 



*V , 'A 



where 6 = 2x,i-f- and c = 2x/-— . In order to find the lepton mixing matrix we need to 
A A 

diagonalize the charged lepton mass matrix and, performing a double expansion in the 
parameters u and t, we get 

mf ag = U} c m e Ui = (y e u 2 t, y^u, y T )uv d , (18) 

where the unitary Ue results to be the unity matrix, apart negligible corrections of order 
u 2 t 2 . As a consequence we get that 

Upmns = U\Utb = Utb ■ (19) 

When we introduce the NLO terms in the Lagrangian and the corrections in the VEVs of 
the flavons, we expect corrections to the TB mixing of relative order u, as we will discuss 
in the next sections. As a consequence, it provides an upper bound on the parameter u, 
indeed the maximum deviation from the TB pattern, which we can accept, is 0.05. For 
u > 0.05 the model provides a 9 12 angle which is not in agreement at 2a error with respect 
to the experimental data. The hierarchy of the charged leptons comes directly from the 
symmetry of the model and it is possible to get a constraint on the parameters u and t: 
indeed, for a very low tan/3 value, the requirement for the Yukawa of the r lepton to be 
in the perturbative regime (y T < 4tt) corresponds to a lower bound for u of about 0.001. 
However, using this value for u, we require a particularly large value for y^ in order to 
fulfill the measured value for the ratio m A1 /m T : from the requirement that also remains 
in the perturbative regime, the lower bound on u is raised and we fix it at 0.01. In order 
to explain the ratio rrie/m^, we get a range of values for the parameter t, which is to be 
similar to that for u. Finally we can write 

0.01 < u,t < 0.05. (20) 



3.2 The Quark Sector 

In this part we illustrate the model in the quark sector, getting a good approximation of 
the experimental quark mixing matrix. In table we can see the quark sector fields of the 
model and their transformation properties under S4 x Z 5 x U(1)fn- The superpotential in 





D q 


<?3 


u c 


d c 




s c 




b c 


e 






A 


f 






2 


u 


h 


h 


h 


1 2 


li 


h 


ii 


3i 


2 


3i 


2 


1 2 


^5 




UJ 3 


1 


1 


UJ 2 


UJ 2 


^ 2 


UJ 2 


i 


UJ 2 


cu 2 


a; 3 




1 


U(1) FN 








2 


1 














-i 


















Tabic 2: Transformation properties of all the fields in the quark sector. 



the quark sector can be written as 



where 



+ Eli |r* c (VP) K + Eli ^ ^ 
+ El! ^ (a^ 2 } ) ^ + YLx ^ {D q xP)'h d + 
+ Eli ^ {D q x?)h u + Eli * c (flrf)' /*+ 

+ El, (Vf >)' />„ + Eli ^fV (D q xP)'h d+ 
+ Eli (^ (4) ) *. + Eli (D q xf) k 

+ Eli ^ + Eli ^ M 5 ^ 



= {7777 + v^} 

^ (2) = {^e', AAA, AA</>, A w , <p<p<p} 

x® = {^a, m e'e'} 

X (4) = { w , AA} 
X*- 5 -* = {-i/^A, ^^(^, rjrjip} . 

With this setting, the mass matrix for the up quarks is 

VuuH 2 y u u 2 t 2 y ut uH 2 
y cu u 3 y c u 2 y ct u 
VtuU 3 VtcU 2 Vt 



and for the down quarks is 



y d ut y d ut y db u 2 t 
y sd u 2 y s u y sb u 
ybdU 2 y hs u y b 



where the Yukawas are the sum of all the different terms, which appear in the superpoten- 
tial. 

These mass matrices can be diagonalized by the following transformations: 



m 



diae 



U ] u cm u U u = (y u u z t 2 , y c u\ y t )v u m^ ag = XJ\ c m d U d = (y d ut, y s u, y b )uv d (24) 



diag 



where the unitary matrices can be written in terms of order of magnitude of u and t as 




u. 



Ua 



1 

-0{u) 
-0{u 2 ) 

1 

-0(t) 
-0(ut) 



0{u) 
1 

-0(u) 
0(t) 

1 

-O(u) 



0(u 2 ) 
0{u) 
1 

0(ut) 
0(u) 
1 



(25) 



The resulting quark mixing matrix is 
/ 



V, 



CKM 



u!u d 



Vs 



Vcu 
Vc 



II 



\ 



VbsVsd 



Vcu Vsd 
Vc Vs 

VbdV, 



VbVs 



II 



a 



1 

ybs 

Vb 



-u 



ybdVc - yb S y c 

VbVc 

Vbs 

— u 

Vb 

1 



u 



(26) 



/ 



In order to fit the experimental values of the mixing angles we need to invoke a moderate 
fine-tuning in some parameters. The (23) entry of Vqkm has to be of order 9 2 ~ 0.05 and 
therefore suggests for u a value close to its upper bound. However this is not a strict 
constraint because this value can be well explained for the entire range of u considering 
the Yukawas. On the other hand, the entry (12) requires an accidental enhancement of 

the combination ( — — ) of order 1 16 C ~ 4 in order to describe the correct Cabibbo 
V. Vs Vc J 

angle. It is possible to explain such an enhancement considering particular values of the 
relative phase, („, between — and which is connected to the CP violating phase: if 

Vs Vc 

( q = it, then the two factors sum up and the required values are easily explained. 



4 The Vacuum Alignment 

In the following we present the mechanism to get the particular VEV alignment used in 
the previous sections. In table [3] we illustrate all the flavon fields of the model and a set 
of new fields, the driving fields, defined as scalar fields with vanishing VEV, which are 
used only to select the particular solutions of the scalar potential. In order to distinguish 
between the matter fields, the flavons and the driving fields we introduce an additional 
U(1)r, under which the fields have quantum number 1, and 2 respectively. The usual 
R-parity, useful to avoid FCNC in the supersymmetric extensions of the SM, is a discrete 
group of this U(1)r. 
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u 4 
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1 


1 



Table 3: Transformation properties of the flavons and the driving fields. 



The driving superpotential is 

w d = g 1 {A A<p)+g 2 {ip AA)+g 3 {<p°ipip)+ 

+h(i>°i>i>) + W°i>ri)+ (27) 

The equations for the minimum of the scalar potential are obtained deriving by the 
driving fields: 

0i(y>iA 2 - p 2 A 3 ) = (28a) 

gi((piA 1 -(p 2 A 2 ) = (28b) 

^(^3-^0 = (28c) 

( ?2 (A 2 + 2A 1 A 2 )+^? = (29a) 

(7 2 (A 2 + 2A 1 A 3 )+^2 = (29b) 

2/i(^ 2 - iMa) + f2M2 + Wfe) = (30a) 

2f x ^l - V^s) + / 2 (77i^i + Wfe) = (30b) 

2/i (^1 - ifak) + hM% + riM = (30c) 

M f (' + /!i( W2 - W )=0 (31) 

The equations can be divided into almost separated groups. The first five equations, 
(I28ap - fl29bj) . are satisfied by the alignment 

(A) = ( 1 j v A (if) = ( } ) v v , (32) 

which is a stable solution of the scalar potential, with 

v\ = —^-v 2 v v undetermined. (33) 

The three equations (I30al) - (j30c|) . almost separated from the others, are satisfied by two 
different patterns: the first is 

W = ( 1 ] H h)=(°i) v v ^ 



with 

v i> = ~~M~ v v v v undetermined (35) 
2/i 

and the second is 




with v v and ty, undetermined. Only the first solution provides the results presented in the 
previous sections and we need of some soft masses in order to discriminate it as the lowest 
minimum of the scalar potential. We manage in doing it, considering some Z 5 -breaking 
soft terms involving if; and 77, which in the most general form can be written as 

m 2 ^\ip\ 2 + m 2 \r]\ 2 + rrv^ipip + fh 2 r}r] . (37) 

Assuming that TT& < the first two terms stabilize the potential for both the vacuum 
configurations. On the other hand the last two terms vanish for the first vacuum configu- 
ration and get a value different from zero in the second one. With an apposite choice of the 
soft parameters, these contributions can be positive, distinguishing the two configurations 
of VEVs and assuring that one in eq.( l34"l) as the setting with the corresponding lowest 
minimum. 

Acting on the configurations of eq. fl32l) or eq. (1341) with elements of the flavour symmetry 
group S4, we can generate other minima of the scalar potential. These new minima are 
physically equivalent to those of the original sets, but it is not restrictive to analyze the 
model by choosing as local minimum exactly those ones in eqs. f)32l) and (|34p (it is possible 
to show that the different scenarios are related by field redefinitions). 

The last equation (13TT) connects all the sectors and fixes the VEV of £' 

(£') = v? = ^v v v v . (38) 

For the flavon field 9, related to the Froggatt-Nielsen symmetry, the non vanishing 
VEV is determined by the D-term associated with the U(1)fn symmetry (see [13] for more 
details). The D-term in the potential is given by: 

V D = \{M 2 FI -g FN \9\ 2 + ...) 2 (39) 

where g F N is the gauge coupling constant of U(1)fn and M FI is the Fayet-Iliopoulos term. 
The vanishing of Vd requires 

g FN \6\ 2 = M 2 FI . (40) 
Assuming that M FI /g FN is positive, this condition fixes the VEV of 9, given in eq. fl20l) . 



5 NLO Corrections 



We now study the deviations to the LO results. We first present the analysis for the VEV 
alignment and then we move to the mass matrices. 



5.1 The VEV Alignment 



Here we only summarize the results for the vacuum alignment, while a detailed study is 
developed in the Appendix B. The part of the superpotential depending on the driving 
fields A , tp°, if; and £'° is modified into 

w d = w° d + Sw d , (41) 

where w® corresponds to eq. ( 127|) and Swd is the most general quartic, ^-invariant polyno- 
mial linear in the driving fields: 

a* = j ( E zj? + E «tf + E *''f + E ) < 42 > 

\i=l i=l i=l i=l / 

where Xj, u>j, s« and Vi are coefficients and j/j A °, If , , if j represents a basis of 

independent quartic invariants (the list of all the J*° are present in Appendix B). The new 
minimum is obtained by searching for the zeros of the F terms, looking for a solution that 
perturbs eq. (j!5|) to first order in the 1/A expansion: denoting the general flavon field with 
$, we can write the new VEVs as 

($ l ) = ($ i )( L °) + ^. (43) 

All the perturbations are non vanishing, a part drji and Sr)2 and one of the perturbations 
in the neutrino sector, which remains undetermined. On the other hand the NLO terms 
fixes the relation between and v v . We can conclude that the VEV alignment in eq. (|15|) 
is stable under the NLO corrections and the deviations are of relative order u with respect 
the LO results. 



5.2 The Mass Matrices 

In this part we present the corrections to the mass matrices due to the higher order terms 
in the matter superpotential and the deviations to the VEV alignment. 

Lepton Sector : the superpotential for the charged leptons can be written as 

wi = w® + 8w e (44) 

where w® corresponds to eq. ffTTT) and 5w£ contains all the NLO terms. We note that 
the LO operators related to e c completely fill in the first line of rri£ and, as a result, 
the corrections can be reabsorbed in the LO parameters. For this reason, we avoid to 
specify the NLO operators of 5we related to e c , reporting only those ones connected 
to /i c and t c : denoting A and <p with $ v and ip and rj with $£, we can write 

C C 

^{t$&& v + mt'?) , ^(t$„$ v $ v + e$t$it') ■ (45) 

These corrections have to be added to those ones originated by w® considering the 
deviations at the NLO to the vacuum alignment. Finally the corrected charged lepton 
mass matrix has the following structure 

0(u 2 t) 0{uH) 0(uH) 
mi = | 0{u 2 ) 0{u) 0(u 2 ) | uv d , (46) 
0{u) 0{u) 0(1) 



where only the order of magnitude of the single entries are reported. As a consequence 
the unitary matrix Ue, which corresponds to the transformation of the charged leptons 
used to diagonalized me, is modified in the following way: 
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T{ 2 u 




T( 2 u 


1 


rpe 
1 23 




-T23U 


1 



(47) 



where the parameters are factors of order one. 



A similar analysis can be performed for the neutrino superpotential 

tu„ = iy[J + 6w v (48) 
where w® corresponds to eq. (fl~2]) and 5w u contains the only NLO operator, 

^(£hJh uV ye . (49) 

In addition to this correction, we have to consider those ones from w®, with the 
deviations at the NLO to the VEVs. As a consequence the neutrino mass matrix is 
corrected by terms of relative order u in every entry. Now the TB pattern has to be 
modified in order to diagonalize m u and we can write 

U v = U TB + SU v u (50) 

where 5U V can be parameterized by three angles, T± 2 , T 2 V 3 and T± 3 , in a similar way 
as in eq. (j4T]) . 



Finally, summarizing all the corrections from the higher order terms, deviations to the 
neutrino mixing matrix of relative order u with respect the LO results are generated. 
The corrected neutrino mixing angles are modified as follows: 

tan£ 23 = -1 - 2u | T 2 e 3 + ^iV^ 27%, ( (gl) 



1 3m 

tan 612 = — 1= - 

V2 4 



(V2 (7T 2 + 2* ) - V3 (Tf 2 + TQ) (52) 
tan 9 13 = ^= (V6 (T& - Tf 3 ) + + 2^ - 3^) . (53) 



We can conclude that the NLO corrections originate deviations to the TB mixing 
angles of order u. 

Quark Sector : the analysis for the up and down quark mass matrices is simpler than 
the previous case, because m u and rrid do not have any vanishing entry at LO and 
therefore the corrections from the NLO operators of the superpotential and from 
the deviations to the VEVs introduce correcting factors of relative order u in each 
entry of the mass matrices. As a result the quark mixing angles receive deviations of 
relative order u, which do not spoil the LO results. 



6 Conclusions 



The aim of a flavour model is getting the correct mixing angles and mass hierarchies of 
both leptons and quarks, without inducing not observed processes, like FCNC and proton 
decays. Moreover in the context of non-Abelian flavour discrete symmetries, we face off 
the further problem of keeping and preserving a different VEV alignment for the flavons 
and also this requirement has to be naturally fulfilled^. All these points are separated one 
from each other and it seems very hard to get all of them at the same time using a single 
flavour symmetry group. However if a model manages in doing it, it will be considered as 
the most promising model in order to describe nature. Trying to understand what is the 
best candidate, many models have been proposed based on a product of different symmetry 
groups. However only few of them appear interesting: we consider fundamental aspects 
the lack, or at least a moderate amount, of fine-tunings, the smallness of the number of 
elements in the complete flavour symmetry group and in the list of the new scalar fields, 
the flavons. 

Following these points, in this paper we have presented a model for fermion masses and 
mixing angles based on the flavour symmetry group S4 x Z§ x U(1)fn- The main aspect is 
the spontaneous breaking of S4, which guarantees the TB pattern as the neutrino mixing 
matrix at LO. This feature is common to other models, like for example those based on the 
group A4 [8-15] or those containing v4 4 , as the group T' [20,21]. The choice of 5*4 has been 
suggested by the recent work by C. S. Lam [22]: S4 results to be the only group (with all 
the groups containing S4) which predicts the TB mixing in a natural way, namely without 
ad hoc assumptions, from the group theory point of view. This result is completely apart 
from all the possible realizations of a model based on £4 and predicting the TB pattern: 
indeed, from the model building point of view, the most economic group which realizes this 
particular neutrino mixing matrix is A4. However, there are other reasons which enforce 
the use of S4: in the ^4-based models, it seems very difficult and unnatural to generate 
the correct mass hierarchies and mixings for quarks. 5*4 represents a viable solution to this 
problem, because it contains a doublet representation more than A4, which can be used in 
order to describe quarks. 

This is not the first attempt in this direction: in [20] the group T' has been studied 
with good results, getting the TB mixing for leptons and a realistic quark mixings together 
to correct mass hierarchies. Unfortunately, this model suffers of a fine-tuning in order to 
generate the up-quark mass and the (12) entry of the CKM matrix: these negative aspects 
have been already underlined in the papers by Barbieri et al, studying the continuous 
group U{2) [19], and it is connected to the fact that the first two quark families, both left- 
and right-handed, transform as doublets. In our model, we followed a different strategy, 
letting only the left-handed quarks of the first two families transform as a doublet, while 
the right-handed transform as singlets of S4. We manage in getting the correct up quark 
mass, but we ask to some parameters to combine in such a way to bring an unjustified 
factor of order 1/9 C or, from an alternative point of view, we ask to the phase ( q to be close 
to 71. From this point of view, we only partially overcome to the problems of the previous 
models based on U (2) and T' . 

In the lepton sector, our model predicts a neutrino mass matrix which can be diag- 
onalized by the TB pattern and a realistic charged lepton mass hierarchy. With respect 
to the ^U-based models, we predict an inverted hierarchy for the Majorana neutrinos and 



6 Some attempts in which the VEV alignment problem in not present can be found in [34]. 



we get some interesting bounds on |m ee |, which dominates the 0^2/3-decay, on the lightest 
neutrino mass and on the sum of the neutrino masses: 

\m ee \ > 14.4 meV \ m u 3 \ > 0.72 meV 

\m v i\ > 89.4 meV . (54) 

i 

Moreover we have a prediction for \m ee \ in function of Am^ (m , r and the phase (, defined 
in section 2, 

-(1 + r)Am 2 atm + yj (Ami) 2 cos 2 C (3(r - l) 2 + (r + l) 2 cos 2 Q sin 2 C . 

Our predictions are quite close to the future experimental sensitivity, which are expected 
to reach the values of 0.090 eV [28] (GERDA), 0.020 eV [29] (Majorana), 0.050 eV [30] 
(SuperNEMO), 0.015 eV [31] (CUORE) and 0.024 eV [32] (EXO). These aspects can dis- 
tinguish our model from all the others using only observables linked to the neutrino oscil- 
lations. Furthermore, other studies, like on some lepton flavour violating precesses, can be 
performed in order to complete this analysis and better characterize our proposal. 
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Appendix A: The Group S4 



The character table of the group 5*4 is 
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Tabic 4: Character table of 1S4. Ci are the conjugacy classes, n the number 
of elements in each class, h the smallest value for which \ — 1- I n the last 
column we have reported an example of the elements for each class. 



The generators, S and T, obey to the following rules 

S 4 = T 3 = (ST 2 ) 2 = 1 (A.l) 
and can be written in the different representations as 
representation li: S — 1, T — 1 

representation I2: S = —1, T = 1 

1 \ m / uj 



representation 2: S — ( ^ ^ 1 , T — , ^ 2 



representation 3i: S 



representation 32= S 
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The 24 elements of the group belong to five conjugacy classes 



Ci 

c 5 



1 

S 2 , TS 2 T 2 , S 2 TS 2 T 2 

T, T 2 , S 2 T, 5 2 T 2 , STST 2 , STS, S 2 TS 2 , S 3 TS 
ST 2 , T 2 S, TST, TSTS 2 , STS 2 , S 2 TS 
S, TST 2 , ST, TS, S 3 , S 3 T 2 . 



In the 2-dimensional representation the elements are 

1 
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uj \ ( uj 2 



uj 2 / ' V 



1 \ / w\ / a; 2 

C4 



4,5 ■ 1 n ' , ,2 



1 / ' V oj 2 ' \ uj 



while for the 3- dimensional representation 3i the elements are 
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and finally for the 3-dimensional representation 32, the matrices representing the elements 
of the group can be found from those just listed for the representation 3i: for Ci 2,3 are the 
same, while for are the opposite. It is connected with the generator S, which changes 
sign in the 3i and 32 representations: the elements in £1,2,3 contain an even number of S, 
while those in C4 5 contain an odd number of it. 



We now report the Clebsch-Gordan coefficients for our basis. In the following we use 
a,i to indicate the elements of the first representation of the product and fa to indicate 
those of the second representation. 



We start with all the multiplication rules which include the 1-dimensional representations: 
li <8> t\ = T) <8> li = t\ with rj any representation 



U® h 
1 2 ®2 

1 2 ® 3i 

1 2 <8> 3 2 
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2 r 



af3 



a(3\ 
-af3 2 



3i 




The multiplication rules with the 2-dimensional representation are the following: 

' li~ ai(3 2 + ot 2 @i 



2® 2 = li 



with 



2 <g> 3! = 3i © 3 2 



with 



1 2 
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a 2 P 2 
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axj3 2 - a 2 (3 3 
OLifiz - a 2 (3i 
a-iPx - a 2 (5 2 

aif3 2 + a 2 f3 3 
otifii + a 2 (3 2 



The multiplication rules with the 3-dimensional representations are the following: 

li ~ ax(3i + a 2 (3 3 + "3/52 



a 2 j3 2 + + a 3 /3i 
0:3(33 + atifo + 012P1 
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Appendix B: The Vacuum Alignment at NLO 

In this section there is the analysis for the corrections to the vacuum alignment introduced 
by the higher dimensional operators. In table [5] there is a summary of the transformation 
properties of the flavons and of the driving fields. 
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Tabic 5: Transformation properties of the flavons and the driving fields. 



The part of the superpotential depending on the driving fields A , ip°, ip° and £'° is 
modified into 

w d = w Q d + 5w d . (B.l) 



The leading order contribution is 

w d = g 1 (A Aip) + g 2 {cp AA) + g 3 (ip (p(p)- 
+/i(^) + / 2 (^V/) + 



and the minimum is 
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v xV = — tyv . 



;b.4) 



I 1 JWf 

The remaining part, Swj, is the most general quartic, ^-invariant polynomial linear in the 
driving fields: 



\i=l i=l i=l i=l J 



(B.5) 



where x i} w i: Sj and Vi are coefficients and |/ A °, If , if , l| | represents a basis of 
independent quartic invariants: 



If = (A°(A^) 3l )'r 
J A ° = (A°(AA) 3l )'r 
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1 1; 
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(B.6) 



if 



r°e'(Av) 



The new minimum for A, 77 and £' is obtained by searching for the zeros of the F 

terms, the first derivative of Wd + Swd, associated to the driving fields A , tp°, ip° and 



We look for a solution that perturbs eq. (lB.3l) to first order in the 1/A expansion: denoting 
the general flavon field with <£>, we can write the new VEVs as 

(* i ) = {$ i )M + 5$ i . (B.7) 

The minimum conditions become equations in the unknown v v and v v . By keeping 
only the first order in the expansion, we see that the equations can be separated into 
different groups: the first five concern only the neutrino sector, the second three only the 
charged lepton one and the last one connects the two sectors. Finally all the perturbations 
are non vanishing, a part Srji and 5r]2 and one of the perturbations in the neutrino sector, 
which remain undetermined. On the other hand the NLO terms fixes the relation between 
v v and v v . We can conclude that the VEV alignment in eq. (lB.3l) is stable under the NLO 
corrections and the deviations are of relative order u with respect to the LO results. 
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